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B ernard Riemann (1826-66) was a student of Carl Friedrich Gauss (1777- 
1855) at the University of Gottingen in Lower Saxony. Since his youth, 
Gauss had conjectured that the distribution of the primes could approximately 
be given by 

N 

ln(iV) ’ 


the approximation becoming better as N became larger. In time, this 
conjecture became known as the prime number conjecture and then as the prime 
number theorem or PNT. Today we write this as 

k(N) ■ 


N 


ln(JV) 


where n(N) represents the number of primes less than a given number N. 
A second way to write the PNT is n(x) ~ Li(jc) where Li(x) is called the log 
integral function and represents the definite integral of ( as t varies from 
2 to x. 

In 1859, on the occasion of being elected as a corresponding member 
of the Berlin Academy, Bernard Riemann presented a lecture with the tide 
On the Number of Primes Less Than a Given Magnitude, in which he presented 
a mathematics formula, derived from complex integration, which gave a 
precise count of the primes on the understanding that one of the terms in the 
formula, which depended on a knowledge of the non-trivial zeros of the zeta 
function, could be evaluated. Riemann had calculated some of the non-trivial 
zeros and found them all to have a real part equal to 0.5. He conjectured 
that every zero of the zeta function had a real part equal to 0.5. This became 
known as the Riemann conjecture which evolved into the Riemann hypothesis as 
more supporting evidence became available. 

Riemann’s 1859 paper can be found in Stephen Hawking’s 2005 
book God Created the Integers. Riemann’s solution is given using complex 
number integration and is assessable to university students studying higher 
mathematics courses. A more assessable solution for school students who are 
studying advanced mathematics (e.g., Specialist Mathematics in the Australian 






















Curriculum: Mathematics) is given injohn Derbyshire’s very readable 2012 book 
Prime Obsession. 

The aim of this paper is to show that the zeros of the Riemann zeta function 


C( s) = V —, where s = a + ib 

^ k=1 k 


all have the real part equal to one half. 

This is demonstrated in two ways or methods. Each method shows that 
there can be only one value for the real part a, and since we know that there 
are zeros with a = ' then the proof is complete. 


Notation 

(a) The Riemann zeta function 


y- |T|=i + ± + l + J- + ...=«,) 
r r 3' r s 


where s = a + ib and a and b are real, 
(b) The alternating Riemann zeta function 

1 


—(—IV 


k-i 


1111 r is 

_ I 7_ F + 3 7 ~4 7 + '" _ ^ 5 


where s = a + ib and a and b are real. 

(c) The symbol ~ between two functions or series means that the two 
are asymptotically equal to one another. One series might represent 
a convergent series in the variable x and the other might represent a 
divergent series in x but the two may be equal to one another to some 
degree, say four decimal places, when a finite number of terms of the 
divergent series is used. An example of this is given below: 

V°° J_ i , Af_ 4 -:®£_/*\ 
~ 2 1 2 2 2 3 2 4 2 5 2 6 '" V ; 


Now this is an example of an asymptotic series. We are approximating a 
slowly converging series with a diverging series. The numbers B 0 , B v B.„ B 3 , 
... represent the Bernoulli numbers 1, — 4> 0, —-ip 0,-L, 0, — JL, 0,..., with 
B 2k+1 = 0for k>l. 

The value of the right-hand side of (*) for varying numbers of non-zero 
terms is (to 6 decimal places): 


One term: 

1 

Two terms: 

1.5 

Three terms: 

1.625 

Four terms: 

1.645833 

Five terms: 

1.644792 

Six terms: 

1.64498 

Seven terms: 

1.644913 

Eight terms: 

1.64495 
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Now the exact value of the left-hand side of (*) is * = 1.644934 (to 6 
decimal places). 

So we see that the right-hand side is equal to the left-hand side correct 
to four decimal places if seven terms of the right-hand side series are used. 
However, if more than 11 terms of the series on the right-hand side are used, 
then the two sides are no longer equal correct to four decimal places. This 
is mainly because the first eight Bernoulli numbers lie between zero and 1.1 
in magnitude while the next eight Bernoulli numbers lie between 1.1 and 
27 298 231.1 in magnitude, which far outpace the denominators of successive 
terms which are only increasing by a factor of 2. 


Method one 

In this method, in s = a + ix, a is considered constant and x is variable and 
C, a (a + ix) is found asymptotically in the form (A 0 + AjX + A 2 x 2 + A 3 x 3 + A 4 x 4 + 
...) + i(C 0 + C 1 x+ C 2 x 2 + C 3 x 3 + C 4 x 4 + ...) where the A ; and C i are real numbers. 
Justification for assuming such a result is as follows: 

C a (5) = (l-2 4 -%s) 

= (l-2 1 - s )X~ - 

= (l-2 1 "') X L ( Y ) [cos (* ■In W) -1 sin (x In (&))] 

Now, cos(xln(&)) and sin(xln(&)) can be written in terms of power series 
in the variable x. 

Also, (1 - 2 1-s ) = (1 - 2 1 ~ (a+!5t) ) can be written as a series in x with real 
coefficients. 

Hence we have an infinite sum of terms of the form .S', (x) + iS 2 (x) multiplied 
by a series in x. This will give an expression in x with real coefficients. This is 
the proposal and justification. 

However, the method to be used is not that oudined above. Instead, we 
equate C, a (s ) and (A 0 + A x x + A 2 x 2 + A 3 x 3 + A 4 x 4 + ...) + i(C 0 + C 1 x + C 2 x 2 + C 3 x 3 
+ C 4 x 4 + ...) with the understanding that they are asymptotically equal. 

So, let = (A) + A x + A* 2 + A^ 3 + A 4 x 4 + ...) 

+ i( C 0 + CjX + C 2 x 2 + C 3 x 3 + C 4 x 4 + ...) (A) 

We wish ultimately to find expressions for the coefficients A ■ and C t . Then we 
can look for solutions to £ 0 (s) = 0 in the domain 0 < x < 1. 

So to begin: 

1111 

-^a+ix ^ a+ix <%a+ix ^a+ix 

= ^Aq + AjX + AqX" + A 3 x 3 + A 4 x 4 +...j + + CjX + C 9 x~ + C 3 x 3 + C 4 x 4 + .. .j (1) 
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Let x = 0 in equation (1) 

-1-h . . . — A {) + iCn 

jfl 2® ^ u 

As the left-hand side is real, the right-hand side is real. So 
C n = 0 and A, =-1- 

U aj -j a ^ a ^ a 

Now use differentiation of equation (1) to find more coefficients. The first 
differentiation gives: 

/(In (2) In (3) In (4) 


I f^a+ix ^ a+ix ' ^ a+ix 

— ^ Aj + 2A^x + 3 A^x^ +... j + i + 2 C 2 X + C%x^ +... j 


Letting x = 0 in equation (2) we obtain 

^ 2 a 3 a 4 a J 11 

Equating real and imaginary parts, we hnd 

„ A ^ (l n (2) ln(3) ln(4) 

A, = 0 and C, =-1-... 

1 JL ^a ^a ^a 

Subsequent differentiation and setting x = 0 gives the following results: 
. , ( 1 1 1 1 1 

0 l 2 “ 3 s 4 a 5° 6“ 

C 0 =0 

Aj = 0 

ln2 ln3 ln4 ln5 
C, =-+-+... 

1 2° 3“ 4“ 5 a 


( 2 ) 


A, = — 

- 2 ! 

C 2 =0 

a 3 = o 


hr 2 In 3 hr 4 hr 5 

-1-h ... 

2 a 2> a 4“ 5“ 


1 

( In 3 2 

In 3 3 

In 3 4 

In 3 5 

3! 

y 2“ 

3“ 

4 a 

5“ 

1 

f In 4 2 

In 4 3 

In 4 4 

In 4 5 

4! 

t 2“ 

3“ 

4 a 

5 a 


+ ... 


+ ... 


C 4 =0 

a 4 = o 


c 5 = — 

5 5! 


In" 2 ln D 3 ln s 4 In 5 5 
- 1-h ... 

2« 4 . a 


A 6 “ 


6 ! 


In 2 hr 3 In 4 In 5 
-+-+ ... 


In general, A n = 0 if n is odd; C n = 0 if n is even; 

n / 1 \ I n 7 

A = (- 1) 2+1 if 71 iS eVel1 alld n * 0 

c - =( - 1 )l "tii) s « < - 1) ‘ 1I ? if ” isodd 
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So the first aim has been achieved, that is, we have found C, a (a + ix) 
asymptotically in the form 

(A 0 + A x x + A 2 x 2 + A 3 x 3 + A 4 x 4 + ...) + i(C Q + C 1 x + C 2 x 2 + C 3 x 3 + C 4 x 4 + ...). 
In fact, noting that some of the coefficients are zero, we can write: 

C, a (a + ix) ~ (A 0 + A 2 x 2 + A 4 x 4 + A 6 x 6 + A g x 8 + ...) + i(C 1 x+ C 3 x 3 + C 5 x 5 + C 7 x 7 + ...) 
Notice that the imaginary part of C, a (a + ix) has x as a factor. 

Now let us go to our second aim which is to consider the equation 
l^ a (a + ix) = 0 for 0 < x < 1. 

So we have 

(A 0 + A 2 x 2 + A 4 x 4 + A 6 x 6 + A 8 x 8 + ...) + t( CjX + C 3 x 3 + C 5 x 5 + C 7 x 7 +...) = 0 + i.O 
Equating real and imaginary parts, we have 

(A 0 + A 2 x 2 + A 4 x 4 + A 6 x 6 + A 8 x 8 + ...) =0 and (C 4 x + C 3 x 3 + C 5 x 5 + C 7 x 7 + ...) = 0. 
Let these two equations be rearranged as follows: 

The first becomes 


A 


f A A v® A ^ 

_ iAi 1 Ao i 1 Ij iAi 1 1/- 

1 + —— +-— +-=2. + ... 


V 


AAA 


= 0 


/ 


and the second becomes 


CjX 


2 . x 2 C 3 x 4 C 5 x 6 C 7 


1 + 

V C i 


Cl 


+ ... 

Cl J 


= 0 


Dividing these last two equations by the non-zero numbers A 0 and C 1 
respectively, we have 


A x n 4 A v* ^ A 

/lo X IAa X /I c 

1 +-- +-- +-2- + ... 

AAA 


r 


= 0 and x 


x 


1 + 

V C i 


2 c v 4 c v 6r 


- + - 


Q c, 


-+... 




Now let us consider solutions of £ a (a + ix) = 0. These must also be solutions 
of the above two equations. Also we know that x = 0 is not a solution of 
C, a (a + ix) = 0. So we conclude that the three equations 

= 0, C, a (a + ix) = 0 


v 2 A A x°^ A 

A 1 At) iA. /I a iA 1 1 A/- 

1 + - ~ + --+ -- + . 

A 

II 

© 

( v 2p „4C v.6 p "\ 

1+ 3 + 5 + 7 +... 

v A A A 

J 


1 c, C, c, J 


have the same solutions and that since the first two series are monic, they are 
equal. 


A _ C3 
A Q 


Equating coefficients of the terms in x, we have 
A C 

Then it follows that —-- = 0 . 

A A 

Note here that it should be remembered that the equal signs being used 
above are not actually equal signs but rather they represent approximations 
to a certain number of decimal places, or asymptotically equal signs. 


A 9 Co 
Now let y = —— — 

A Q 

and using the values of the coefficients found previously, we have 


f In 2 2 

In 2 3 In 2 4 

In 2 5 

A 


( In 3 2 

In 3 3 In 3 4 

In 3 5 ^ 

{ 2 “ 

3 ° ‘ 4 “ 

5 “ + ' 

J 


[ 2 “ 

1 

5 “ + "'J 


y UuAA + i L + i_L + ^ + UJ fln2_ln3 + ln4_ln5 + i 

A _ 2^ + 3“ _ 4 7+ fA~6“ + "’J [Y a ly r + ^i a 5 r + "'J 
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Notice that y has been expressed in terms of one unknown, a. This can 
be shown on a graph using the software Graphmatica, but first the unknown 
must be replaced by the variable x. This graph is shown below for the domain 
0 < x < 1, that is, 0 <a< 1. Each of the four series involved has been truncated 
to ten terms. 

The graph shown in Figure 1 is as it would appear on a hard copy using the 
Graphmatica software (see Appendix 1 for a larger section of this graph). 



Figure 1. 


Equations on screen: 

y=.5*(f2(x))/(fl(x))+l/6*(h2(x))/(hl(x)) 


Functions used by these equations: 
f=l 

fl(x)=(2 A -x)*(ln(2)) A (2*f-2)-(3 A -x)*(ln(3)) A (2*f-2) 

+(4 A -x)*(ln(4)) A (2*f-2)-(5 A -x)*(ln(5)) A (2*f-2)+(6 A -x)*(ln(6)) A (2*f-2) 
-(7 A -x)*(ln(7)) A (2*f-2)+(8 A -x)*(ln(8)) A (2*f-2)-(9 A -x)*(ln(9)) A (2*f-2) 
+(10 A -x)*(ln(10)) A (2*f-2)-(ll A -x)*(ln(ll)) A (2*f-2) 
f2(x)=(2 A -x)*(ln(2)) A (2*f)-(3 A -x)*(ln(3)) A (2*f)+(4 A -x)*(ln(4)) A (2*f) 
-(5 A -x) * (In (5)) A (2*f)+(6 A -x) * (In (6)) A (2*f) 
-(7 A -x)*(ln(7)) A (2*f)+(8 A -x)*(ln(8)) A (2*f) 

-(9 A -x)*(ln(9)) A (2*f)+(10 A -x)*(ln(10)) A (2*f)-(ll A -x)*(ln(ll)) A (2*f) 
hi (x)=(2 A -x) * (In(2)) A (2*f-l)-(3 A -x) * (In (3)) A (2*f-l) 

+(4 A -x)*(ln(4)) A (2*f-l)-(5 A -x)*(ln(5)) A (2*f-l)+(6 A -x)*(ln(6)) A (2*f-l) 
-(7 A -x) * (In(7)) A (2*f-l)+(8 A -x) * (In (8)) A (2*f-l) 
-(9 A -x)*(ln(9)) A (2*f-l)+(10 A -x)*(ln(10)) A (2*f-l)-(ll A -x)*(ln(ll)) A (2*f-l) 
h2(x)=(2 A -x)*(ln(2)) A (2*f+l)-(3 A -x)*(ln(3)) A (2*f+l) 
+(4 A -x)*(ln(4)) A (2*f+l)-(5 A -x)*(ln(5)) A (2*f+1) 
+(6 A -x)*(ln(6)) A (2*f+l)-(7 A -x)*(ln(7)) A (2*f+1) 
+(8 A -x)*(ln(8)) A (2*f+l)-(9 A -x)*(ln(9)) A (2*f+1) 

+(10 A -x)*(ln(10)) A (2*f+l)-(ll A -x)*(ln(ll)) A (2*f+l) 
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In the following conclusion, use is made of the fact that £(s) and £ a (s) 
have the same zeros, that is, £(s) = 0 and £ a (s) = 0 have the same roots, e.g., 
s = 0.5 + il4.134725 approximately. 

Conclusion for Meth od 1 

From the graph we see that there is only one intercept in the domain 0 < x < 1; 
this occurs at about x = 0.45. This means that there is only one possible value 
of a in the domain 0 < a < 1 at about a = 0.45. Now it is known that all non¬ 
trivial Riemann zeta zeros have real part a lying in the interval 0 < a < 1. But 
our graph shows that there is only one possible value for a in this interval. And 
since it is also known that there are Riemann zeta zeros with real part a = 0.5, 
then we conclude that the value of a = 0.45, shown on the above graph, is 
an approximation for a = 0.5. Then we finally conclude that if s = a + ib is a 
solution to the equation £(s) = 0 or t, a (s) = 0 then a = 0.5. 


Method 2 


The second method begins with the asymptotic formula 

y- f n \ , jo By , B 9 _s , B 4 s(s + 1)(s + 2) 

^ k =\k s )~ 2 i _ 1 ( 5 - 1 ) 2' 5 2 ! 2 S+1 4 ! 2 f+3 

, B 6 s(s + 1)(s + 2)(s + S)(s + 4i) , 

eiiF (B) 

where B { are the Bernoulli numbers and s = a + ib. 

In this method, the intention is firstly to truncate the right-hand side of 
the above equation to two terms, let s = a + ix as in Method 1, express the 
truncated series in the form S 1 (x) + iS 2 (x), where 

Sj (x) = (D 0 + D 1 x + Z) 2 x 2 + D 3 x 3 + D 4 x 4 + ... ) , 

S 2 (x) = (E 0 + EjX + E 2 x 2 + E 3 x s + £ 4 x 4 + ...), 
compare the coefficients of x 2 as in Method 1 and produce a graph. Secondly, 
truncate the right-hand side of (B) to three terms and do the same as above 
to produce a graph; then truncate the right-hand side to four terms, and then 
five terms, producing four graphs in total. Finally, a conclusion is arrived at in 
a similar vein to that in Method 1. 


Step 1 

Truncate the right-hand side of (B) to two terms and let s = a + ix. 




■ 1 + 


B n 




= 1 + 


B 0 2 1 

s-l 


= 1 + 


BgT 


a + ix -1 


1 + B 0 {e In2 ) 1_S - {a 1MX - 

v 7 (a—l + ix) (a—l — ix) 

1 + Bg(e™ f~ aHX -- 

v 7 (a-\ + ix) (a-l-ix) 
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Now multiply both sides by [ (a - 1)~ + x 2 ] to obtain 

[(a — 1)'’ + —'j ~ [(a —1)" +x“] + 5 0 2 1 " [(a-l)-tx][cos(xln2)-tsin(xln2)] 


Writing 

and 


cos(xln2) = 1 — 


x l In 2 2 x 4 In 4 2 


2 ! 


4! 




. , . _ x 3 In 3 2 x 5 In 5 2 

sin (x In 2) = xln2-—-t-- 


3! 


5! 


we arrive at the following: 
[(a-l) 2 +x 2 ] C, a (a + ix) 


-|[(«-1) 2 + B 0 2 1 ->-1)] + l 
+ i\ -Bo 2 1 ”" [1 + (a- 1)In2]x + B 0 2 1- ' 


In 2 2 

l-5 0 2 1_a (a-l)- --B 0 2^ a ln2 x 2 + terms in x 4 , x 6 , ...] 


2 ! 

In 2 2 (a —l)ln 3 2 


2 ! 


3! 


3 -57 

x' + terms m x , x , 


We note that the imaginary part of the above has a factor of x. When this is 
factorised, the result can be written in the form: 

[(a-1) 2 + x 2 ]( l a (a + ix) ~ (D 0 + Z) 2 x 2 + Z) 4 x 4 + .. ,)+ix(^E 0 + f? 2 x 2 + U 4 x 4 +...) 


and then as 

[(ffl-1) 2 + x 2 ]i^ a (« + ix) ~ D 0 


^ Z) 2 x 2 D.x 4 

1 + — 1 — + —— + ... 


A A 


+ ixEn 


f 2 t -1 4 ^ 

, U 9 X iux 


/ 


1 + 


A A 




Using similar reasoning to that used in Method 1, we can say the following: 
The equation 

[(a-1) 2 + x 2 ]i^ a (a + ix) = 0 


implies that 


and 


That 

[(a-1) 2 + x 2 ](^ a (a + ix), 

have the same zeros; 

That 

( 


< A* 2 D 4 x 4 A 

1 + —*— + —— + ... 

A A 


1 ! A x ! A x ! 
A A 


, f) 9 x 2 Z) 4 x 4 

1 + —-+ —-+ .. 

A A 


= o 


= o 


and 


T - ' 2 T - ’ 4 

, h 9 x T 4 x 

1 + —— + —— + .. 

A A 


1 + Af + Ar + - 

y D 0 x D 0 x j 


Equating the coefficients of x 2 in the two series above gives 


y t Q X t Q X j 


and consequently 


Dc, _ E 2 

A A 


Al _ Al _ o 
A A 
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As in Method 1, please note here that it should be remembered that 
the equal signs being used above are not actually equal signs but rather 
they represent approximations to a certain number of decimal places, or 
asymptotically equal signs. 


Let y = 


Do 


D n 



Writing this out fully, we have: 

In 2 2 (a-l)ln 3 2 
2! 3! 

[(«-l) 2 + £ 0 2 1 ^(«-1)] ^ -£ 0 2 w '[l + (a-l)ln2] 


l-5 0 2 1 ^(«-l) h ^ 2 -fi 0 2 1 "ln2 


al-a 


This represents an equation in one unknown, a. 

As mentioned previously, Graphmatica requires us to write equations in terms 
of x. 

The graph of this equation using Graphmatica is shown in Figure 2 as the 
graph y = hi (x) - h2(x). 



Figure 2 


Equations shown in Figure 2: 

y = hi(x) - h2(x) 

Functions used by these equations: 
fl (x)=(x-l) A 2+bO*(x-l)*2 A (l-x) 

f2(x)=l-bO*(2 A (l-x)*(x-l)*l/2*(ln(2)) A 2-bO*ln(2)*2 A (l-x) 
hi (x)=-f2(x)/fl (x) 

f3(x)=l*bO*(2 A (l-x))*(l+(x-l)*ln(2)) 

f4(x)=b0*(2 A (l-x))*(l/2*(ln(2)) A 2+(x-l)*l/6*(ln(2)) A 3) 

h2(x)=f4(x)/f3(x) 
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From the graph we see that there is only one intercept in the domain 
0 < x < 1; this occurs at about x = 0.44. That is, there is only one value of a in 
the domain 0 < a < 1. As concluded in Method 1, this value of a = 0.44 is an 
approximation for a = 0.5. 

Step 2 

Truncate the right-hand side of (B) to three terms. 

The method used in Step 1 above now results in the following equation and 
graph: 

_ /2(a) + /6(a) /4(a) + /8(a) 
y ~ /1(a) + /5(a) /3(a) + /7(a) 

where 

/1(a) = (a-1) 2 + b 0 2 1 ~ a (a-l) 

/2(a) = 1 - b 0 2 1 - 1 ' (a -1) In 2 2 + In 2 

/3(a) = -6 0 2 1_ “(l + (a-l)ln2) 

/4(a) = b 0 2 1 -* [^ln 2 2 + ^(a - l)ln 3 2 
/5(a) = -&j2 _a (a-l) 2 
/6(a) = -^2-^l-^(a-l) 2 ln 2 2j 
/7(a) = -b 1 2~ a (a-1) 2 In 2 
/8(a) = -b x 2”" ^(a -1) 2 In 3 2 - ln21 

From Figure 3, whose equation is y = h7(x) - h8(x), we see that there are 
two values of a as intercepts. One at a = -0.1 approximately and the other at 
about a = -1.3. The important point to notice for this paper is that there is no 
value of a in the domain 0 < a < 1. 
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Step 3 

Truncate the right-hand side of (B) to four terms. 

This results in the graph shown in Figure 4 with the equation y = h 13(x) — hi4(x ). 
Here we find only one intercept at about a = -0.07. There is no value for a in 
0 < a < 1. 



Step 4 

Truncate the right-hand side of (B) to 5 terms and obtain the graph 
y = hl9(x) - h20(x) as shown in Figure 5. Here we find two intercepts at 
a = -0.07 and a = -3.6 but no values for a in 0 < a < 1. 



See Appendix 2 for the full list of equations used by Graphmatica in the 
graphs above. 
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Conclusion for Method 2 

From the four graphs above we see that there is at most one intercept in the 
domain 0 < a < 1; this occurs at about a = 0.44. Now it is known that all non¬ 
trivial Riemann zeta zeros have real part a lying in the interval 0 < a < 1. Our 
graphs show that there is only one possible value for a in this interval. And 
since it is also known that there are Riemann zeta zeros with real part a = 0.5, 
then we conclude that the value of a = 0.44 is an approximation for a = 0.5. 
Then we finally conclude that if s = a + ib is a solution to the equation £(s) =0 
then a= 0.5. 


Overall conclusion 

It has been shown that if ^(a + ix) =0 then there can be only one value of a 
and the value of a has been shown to be of the order of 0.45. It has also been 
shown that since it is known that s = 0.5 + il4.134725... is one solution of the 
equation £(s) = 0 then the real part of any Riemann zeta zero is a = 0.5. 


Appendix 1 


The following is a larger section of the graph used in Method 1 whose equation 
is: 


7 = 


fin 2 2 In 2 3 In 2 4 In 2 5 

-1-r 



^Inf2_lnf3 In 3 4 In 3 5 
1 [ 2“ " 3“ 4" " 5“ 

3 ! fln2 ln3 In 4 ln~5 T 

l~2“ 5^ + "'J 


and its Graphmatica equation is y=.5*(f2(x))/ (fl (x))+l/6*(h2(x))/(hi (x)). 



y=.5*(f2(x))/(fl(x))+l/6*(h2(x))/(hl(x)) 
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Appendix 2 


Graphmatica equations used by the graphs in Method 2. 


bO=l 

bl=-l/2 

b2=l/6 

b4=-l/30 

fl (x)=(x—1) A 2+b0*(x-1)*2 A (1-x) 

f2(x)=l-b0*(2 A (l-x))*(x-l)*l/2*(ln(2)) A 2-b()*ln(2)*2 A (l-x) 
hi (x)=f2(x)/fl (x) 

f3(x)=-l*bO*(2 A (l-x))*(l+(x-l)*ln(2)) 
f4(x)=b0*(2 A (l-x))*(l/2*(ln(2)) A 2+(x-l)*l/6*(ln(2)) A 3) 
h2(x)=f4(x)/f3(x) 
f5 (x) =-l *bl*(2 A (-x)) * (x-1) A 2 
f6(x)=-l*bl*(2 A (-x))*(l-(x-l) A 2*.5*(ln(2)) A 2) 
f7 (x) =-l *bl*(2 A (-x)) *ln (2) * (x-1) A 2 
f8(x)=-l*bl*(2 A (-x))*(l/6*(x-l) A 2*(ln(2)) A 3-ln(2)) 
h3 (x) =fl (x) +f5 (x) 
h4(x)=f2(x)+f6(x) 
h5 (x) =f3 (x) +f7 (x) 
h6(x)=f4(x)+f8(x) 
h7(x)=h4(x)/h3(x) 
h8(x)=h6(x)/h5(x) 
f9(x)=b2*l/2*(2 A (-x-l))*x(x-l) A 2 
fl0(x)=b2*l/2*(2 A (-x-l))*(-l*x*(x-l) A 2*(ln(2)) A 2*l/2 
+ln(2)*(x-l) A 2+x) 

fll(x)=b2*l/2*(2 A (-x-l))*(-l*x* (x-l) A 2*ln(2)+(x-l) A 2) 
fl2(x)=b2*l/2*(2 A (-x-l))*(l/6*x*(x-l) A 2*(ln(2)) A 3 
-l/2*(x-l) A 2*(ln(2) ) A 2-x*ln(2)+l) 
h9 (x) =h3 (x) +f9 (x) 
hi 0 (x) =h4 (x) +fl 0 (x) 
hi 1 (x) =h5 (x) +fl 1 (x) 
hi 2 (x) =h6 (x) +fl 2 (x) 
hi 3 (x) =hl 0 (x) /h9 (x) 
hl4(x) =hl 2 (x) /hi 1 (x) 

fl3(x)=b4*l/24* (2 A (-x-3)) * (x-1) A 2* (x A 3+3*x A 2+2*x) 
fl4(x)=b4*l/24*(2 A (-x-3))*(-3*(x-l) A 2*(x+l) 

-1/2* ((ln(2)) A 2) * (x A 3+3*x A 2+2*x) * (x-1) A 2 
+ln (2) * (3*x A 2+6*x+2) * (x-1) A 2+ (x A 3+3*x A 2+2*x)) 
fl5(x)=b4*l/24* (2 A (-x-3)) * ((3*x A 2+6*x+2) * (x-1) A 2 
-In (2) * (x A 3+3*x A 2+2*x) * (x-1) A 2) 
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fl6(x)=b4*l/24*(2 A (-x-3))*(-l*(x-l) A 2+3*ln(2)*(x+l)*(x-l) A 2 
+((x-1) A 2) *1/6* (x A 3+3*x A 2+2*x) * (ln(2)) A 3 
-((In (2)) A 2) * 1 /2* (3*x A 2+6*x+2) * (x-1) A 2+ (3*x A 2+6*x+2) 

-ln(2)*(x A 3+3*x A 2+2*x)) 
hi 5 (x) =h9 (x) +fl 3 (x) 
hl6(x)=hl0(x)+fl4(x) 
hl7(x)=hll (x)+f!5(x) 
hi 8 (x) =hl 2 (x) +fl6 (x) 
hl9(x)=hl6(x)/hl5(x) 
h20(x)=hl8(x)/hl7(x) 
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